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Abstract: We perform an integral reduction for the 3-loop effective gauge coupling and 
screening mass of QCD at high temperatures, defined as matching coefficients appearing 
in the dimensionally reduced effective field theory (EQCD). Expressing both parameters in 
terms of a set master (sum-) integrals, we show explicit gauge parameter independence. 
The lack of suitable methods for solving the comparatively large number of master integrals 
forbids the complete evaluation at the moment. Taking one generic class of masters as an 
example, we highlight the calculational techniques involved. The full result would allow to 
improve on one of the classic probes for the convergence of the weak-coupling expansion at 
high temperatures, namely the comparison of full and effective theory determinations of the 
spatial string tension. Furthermore, the full result would also allow to determine one new 
contribution of order 0[g 7 ) to the pressure of hot QCD. 
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1 Introduction 

Thermal QCD at high temperatures (T) exhibits three different momentum scales. It has 
been known [1, 2] for a long time that the "soft" static color-electric modes p ~ gT, where g 
is the gauge coupling, are responsible for the slow convergence whereas the "ultra-soft" static 
color-magnetic modes p ~ g 2 T cause the well-known perturbative breakdown [3]. However, 
perturbation theory restricted to the "hard" scale p ~ 2irT can be treated with conventional 
weak-coupling methods, while the soft and ultra-soft scales are only accessible through im- 
proved analytic methods or non-perturbatively via lattice simulations, as is especially the 
case for the ultra-soft g 2 T scale. Here p denotes the characteristic momentum scale, g the 
gauge coupling and T the temperature. The infrared problems which cause the breakdown 
of perturbation theory can be isolated into a three-dimensional (3D) effective field theory 
called magnetostatic QCD (MQCD) and studied non-perturbatively with lattice simulations. 
Before computing various quantities in this framework a number of perturbative "matching" 
computations are necessary [4, 5], in order to relate the parameters of the effective theory 
with those of thermal QCD. 
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The plan of this paper is the following. In Section 2 we review the most important facts 
of the dimensionally reduced effective field theory framework and show how to systematically 
determine the effective gauge coupling gs and screening mass mE- In Section 3 we explain 
some technical details about the integral reduction step, while in Section 4 we discuss the 
structure of the explicit result for the one-, two-, and three- loop corrections, whose rather 
lengthy coefficients are detailed in the Appendix. Section 5 contains the evaluation of a new 
class of master sum-integrals that appear in our result. We finally discuss possible applications 
of our results in Sec. 6, before we conclude in Section 7. 



2 Effective gauge coupling and screening mass 

We consider QCD at finite temperature with the gauge group SU(A r c ) and iVf massless flavors 
of quarks. Before gauge fixing, the bare Euclidean Lagrangian in dimensional regularization 
reads 

Sqcd = J ' dr J d d x£ QC D , (2.1) 

Cqcd = \f; u F; v + fy M Z>„V> (2-2) 

where T is the temperature; d = 3 — 2e denotes the number of spatial dimensions, such that 
Greek indices run as v = 0, . . . , d; F*, = d^-d^A^+gf^A^Al and = ld fl -igA^T a , 
where the T a are hermitian generators of SU(iV c ) with normalization Tr[T a T 6 ] = 5 ab /2; we 
use hermitian Dirac matrices 7/1 = 7^, {~l^,1v} = 25^; g is the bare gauge coupling; and 
tp carries Dirac, color, and flavor indices. For the group theory factors, we us the standard 
symbols C A = iV c , C F = (N 2 - l)/(2iV c ). 

At sufficiently high temperatures, the long-distance physics of Eq. (2.2) can be described 
by a simpler, dimensionally reduced effective field theory [1, 2, 4, 5]: 

<Seqcd = J d d x£EQCD , (2.3) 
£eqcd = \f*F* + Tr[A, B } 2 + m| Tr[B 2 ] + A^Tr^ 2 ] 2 + Ag } Tr[i? 4 ] + . . . , (2.4) 

where i = 1, . . . , d, F§ = dj£?" - + g E f abc B^B^ and A = $ - ^eA- The electrostatic 
gauge fields Bq and magnetostatic gauge fields Bf appearing in the theory above can be 
related (up to normalization) to the zero modes of At of thermal QCD in Eq. (2.2). 

The effective parameters in Eq. (2.4), which we are ultimately interested in, can be 
obtained by matching. This means, we require the same result on the QCD and EQCD side 
within the domain of validity. A convenient way to perform the matching computation is 
to use a strict perturbation expansion in g 2 . On both sides, the expansion is afflicted with 
infrared divergences. These divergences are screened by plasma effects and can be taken 
into account (at least for electrostatic gluons) by resumming an infinite set of diagrams. 
Screening of magnetostatic gluons is a completely non-perturbative effect. For the matching 
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computation, it is not necessary to worry about the infrared divergences because the matching 
parameters are only sensitive to the effects of large momenta. All infrared divergences which 
occur can be removed by choosing a convenient infrared cutoff. It is essential to choose the 
same infrared cutoff in both theories. 

2.1 Relation for m E 

In order to establish a relation between the parameters of the theories Eqs. (2. 2), (2. 4), consider 
the electric screening mass m c \, defined in the full theory 1 by the pole of the static Aq 
propagator, 

= p 2 + U 00 ( P % Q=Qp2= _ m2ci . (2.5) 

On the effective theory side, the electric screening mass is, equivalently, defined as the pole 
of the 3d adjoint scalar Bq propagator, 

= p 2 + m 2 + n EQCD (p 2 )| a , (2.6) 

where IIeqcd denotes the Bq self-energy on EQCD side. 

Noting that the self-energies start at one-loop order, the leading-order solutions for m 2 j 
will be suppressed by the respective coupling parameters, such that p 2 is to be regarded 
perturbatively small, hence allowing for a Taylor expansion of the "on-shell" self-energies 
around zero. For Eq. (2.5), one needs (let us write He = IIoo from now on) 

U E {-m 2 cl ) = He(0) - m 2 cl U' E {0) + ... 

oo oo 

= £ 9 2n n E n(0) - m 2 £ g 2n n' En (0) + . . . , (2.7) 

n=l n=l 

where in a second step we have introduced the n-loop self-energy coefficients II En . From 
Eqs. (2.5) and (2.7), we can express the electric screening mass in terms of Taylor coeffi- 
cients up to next-to-next to leading order (NNLO) 

m^ = 5 2 n E1 (o)+/ [n E2 (o)-n E1 (o)n E i(o)] + 5 6 [n E3 (o) -n E1 (o)n E2 (o) - 

-n E2 (o)n E1 (o) + n^(o)(n E1 (o)) 2 + n E1 (o)(n E1 (o)) 2 ] +o{ g 8 ) . (2.8) 

Diagrams contributing to the various orders of II are depicted in Fig. 1. 

To complete the matching computation for m 2 v we have to compute II E qcd on the 
EQCD side in a strict perturbative expansion. Again treating the "on-shell" momentum p 2 
(as well as the tree-level mass m E ) as perturbatively small, due to the fact that the only scale 
in II E qcd(p 2 ) is p 2 , after Taylor expansion the dimensionally regularized integrals (being 
scale- free) vanish identically 2 . From Eq. (2.6) it hence follows that 

m\ = m 2 cl . (2.9) 

1 In the presence of an infrared cut-off; otherwise, a non-perturbative definition is needed. 
2 Note that this is not the case for the coefficients of Eq. (2.7), since those are vacuum sum-integrals in the 
full theory and hence know about the temperature scale T. 
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Figure 1. The 1-loop, 2-loop and some 3- loop self-energy diagrams in the background field gauge. 
Wavy lines represent gauge fields, dotted lines ghosts, and solid lines fermions. 

2.2 Relation for g\ 

In order to relate the effective 3d gauge coupling g 2 , to the parameters of the full theory, we 
can choose whether to go through a 3-point or a 4-point function, in addition to a 2-point 
function. However, it is further possible to simplify this task to a single 2-point calculation 
using the background field gauge method (see e.g. Ref. [6]). Let us give the main argument 
here, closely following Ref. [7]. 

The effective Lagrangian Eq. (2.4) follows from integrating out the hard (p ~ T) scales 
which, symbolically, produces an expression of the form 

£ cff ~ c 2 {dBf + c 3 g(dB)B 2 + c i9 2 B A + . . . , (2.10) 

where B denotes the background field potential and the coefficients c, = l+0(g 2 ) . Redefining 
now the effective field as B 2 S = C2B 2 , from C c s ~ (dB c g) 2 +czc 2 i ^ 2 ^(dB^B^+c^ 2 g 2 B A S + 
... we can read off the effective gauge coupling (considering the gauge invariant structure 
F 2 ) g c ff = C3C 2 3 ^ 2 g = cl^c^g. Furthermore, since the effective action is gauge invariant 
with respect to both B e ^ as well as B [6], we have 02 = 03 = 04. Finally transforming to 3d 
notation, scaling the fields B — ^T 1 / 2 B 2 and comparing f^ T dr £qcd with £eqcd, it follows 
that 



fjE 



T^c'^g. (2.11) 



Now we proceed in the same way with the effective gauge coupling <?e as for the screening 
mass mE- From Eq. (2.11) we thus obtain 
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9e 



T \g 2 - 5 4 n T1 (0) + g & \{W T1 (0)) 2 - 11^(0) 



+ 



+5 C 



2n^(o)n^ 2 (o) - (n^(o)) 3 - n^ 3 (o) 



+ o 



(a 10 ) } , 



where IIt denotes the transverse part of the (spatial part of the) self-energy 



Ily(p) 



PiPj 
P 2 



n T ( P 2 ) + ^n L (p 2 ) 
p 2 



(2.12) 



(2.13) 



To understand the split-up of II^ in more detail, note that we can choose the external 
momentum p purely spatial, p = (0, p), while the rest frame of the heat bath is time-like, 
with Euclidean four-velocity u = (1, 0), such that u-u = 1, u-p = 0. In this case H^ u has three 
independent components (Iloi, Ho vanish identically). The loop corrections to the spatially 
longitudinal part I1l also vanish (which we will however explicitly check in our computations) , 
such that only two non-trivial functions, ITe and IIt, remain (recall IIe = IIoo)- 

Noting that the class of background field gauges still allows for a general gauge parameter 
£ (we denote (£) he re = 1 - (f Standard), we use the gauge field propagator 



D±(q) 



sab 



(,2)2 



(2.14) 



and verify gauge parameter cancellation in the end of our computations. 



3 The reduction 



After the Taylor expansion and decoupling of scalar products with external momentum, all in- 
tegrals that contribute to the self-energies up to three-loop order that are needed for Eqs. (2.8) 
and (2.12) can be written as 



a,6,c,d,e,/;ci,C2,c 3 



(WTO (P 3 )c 



w 3 \pf\ a \pi\ h \piY m - P2?\ d i(Pi - p 3 ) 2 ] e m - PzW 



, (3.1) 



where Pf = (P 4 )q + p 2 = [(2rtj + c^ttT] 2 + p 2 for i G {1,2,3} are bosonic (fermionic) loop 
momenta for a = (1). The sum-integral symbol in Eq. (3.1) is a shorthand for 



Po 



d d p 
(2^p' 



(3.2) 



where [i is the minimal subtraction (MS) scheme scale parameter, and we take d = 3 — 2e. 

An essential part of this work deals with the reduction of integrals of the type in Eq. (3.1) 
to a small set of master integrals. We use the well-known integration by parts (IBP) identi- 
ties and identities following from exchanges of integration variables. Both are implemented 
in a Laporta algorithm [8] using FORM [9]. Compared to the well-established Laporta-type 
algorithms for zero-temperature reductions, one of the main differences here is that the IBP 
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Figure 2. Top row: non-trivial vacuum topologies at 1-loop, 2-loop and 3-loop. Bottom row: types 
of bosonic and fermionic master integrals. Lines (arrow-lines) corresponds to bosonic (fcrmionic) 
propagators, respectively. 

relations act only within the continuum (spatial) part of our sum- integrals. Another impor- 
tant difference is that in general, linear shifts or exchanges of integration momenta can cause 
a flip of bosonic and fermionic signature of the loop momenta, such that extra care must be 
taken for topology mapping. A precursor of this reduction algorithm had already been tested 
in Ref. [7]. 

The main difference between the outcome of the 1-loop and 2-loop calculation on the 
one side and the 3-loop correction on the other side is that the former ones are expressible 
in terms of 1-loop tadpole sum-integrals which are known explicitly, see App. A. This is no 
longer the case at 3-loop order. The mercedes- and spectacles topology shown on the first 
line of Fig. 2 can be expressed in terms of basketball-type sum-integrals as well as products 
of 1-loop tadpoles. 

4 Structure of the result 

After reduction, we can express all quantities as a sum of 1- and 3-loop master integrals (there 
are no master integrals at 2-loop order, see [5, 10]) of the generic types depicted on the second 
row of Fig. 2, the structure being 

U 3 = ^2a i A i + J2bjB j , (4.1) 

i 3 

where A { = I ■ I ■ I with / G } (4.2) 

and Bj = basketball G {J,K,L} . (4.3) 

A detailed version is given in the Appendix, cf. Eqs. (C.14) and (C.15). 

We have performed a number of cross-checks to confirm the validity of our results: the 
longitudinal parts of the self-energy vanish identically 

n L3 = n' L3 = o for (4.4) 

and the specific combinations of (bare) self-energy coefficients that build up (cf. Eq. (2.8)) 
and g\ (cf. Eq. (2.12)) are gauge-parameter independent up to three-loop order. 
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The one-, and two-loop calculations have already been performed in Ref. [7] which we 
use as another serious cross-check of our independent calculation. We obtain full agreement 
when comparing our Eqs. (C.1)-(C.4) and Eqs. (C.7)-(C.10) with that reference. 

There is considerable experience of how to calculate the genuine 3-loop integrals Bj up 
to the constant term (which can typically only be represented in terms of two-dimensional 
parameter integrals and evaluated numerically), see [11-13]. In Section 5, we add to this 
available knowledge a specific class of 3-loop (basketball-type) sum-integrals which appear in 
our reduced expressions Eqs. (C.14) and (C.15). 

It turns out, however, that most of the pre-factors bj are singular when expanded around 
d = 3 — 2e dimensions. Hence, we need to expand the integrals Bj beyond their constant 
term (in fact, to 0(e) for He and to 0(e 2 ) for iTp). As the conventional techniques for 
computing these basketball-type integrals rely on a careful subtraction of sub-divergences on 
a case-by-case basis, it appears quite difficult to extend the known techniques in order to 
evaluate higher terms in the epsilon expansion. 

To make progress, it might be advantageous to perform a change of basis, see e.g. [14], 
in order to avoid or at least reduce the number of divergent pre-factors. Due to the large 
number of integrals contained in our reduction tables, an algorithmic approach trying out all 
possible different combinations of basis elements might be somewhat involved, but certainly 
possible. 

5 Evaluation of classes of master sum-integrals 

After the successful reduction step, a number of non-trivial three-loop master sum-integrals 
will have to be evaluated. Noting that all bosonic and fermionic one-loop sum-integrals 1^ and 
that appear in Eqs. (C.l)-(C.IO) as well as in Eqs. (C.14), (C.15) are known analytically 
(see App. A), and noting that furthermore all 2-loop structures have been reduced to products 
of 1-loop integrals, let us tackle the first non-trivial sub-class of master integrals, the bosonic 
basketball 

R _ T M,o,o _ vf Qff_ 

&n,m - ^,1,0,0,1,1; 0,0,0 " f J pQR [QY (p _ Q) 2 R 2 (p _ R) 2 > V> A ) 

with N,M > 2. After a careful subtraction of all UV and IR divergences (for more details 
see [11, 12, 15, 16]) we can write Eq. (5.1) as 

B N , M = (3 [A(N, e, 1)S M , + + A° 4+ £ ] + B% M + 

+ 21? S(N, 1, 1; M, 0) + ^ [Q2^p^Q)2 6 m,o + B%8 Mfi + < M , (5-2) 

where /3 = G(l, 1, d + 1) stands for the 4d massless 1-loop bubble and j3 = G(^p, 1, d + 1) is 
the 4d 1-loop propagator, where the function G reads (si2 = s\ + S2 etc.) 

„, «-,>,»,-*/_!_ T(% - ai )T(i - s 2 )T( ai2 - j) 
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and S stands for the two-loop tadpole at finite-temperature, 

S( S1 , s 2 , s 3 ; ai ,a 2 ) = ^ [p2]si [Q2]s2 [{p _ Q)2]S3 

= Is 12S - ai2 /2-i e i( s i> s 2,S3,a 1 ,a 2 ,d) , (5.4) 

i 

where the coefficients ej follow from IBP relations (for an example, see Eq. (5.18) below). 
Furthermore, the abbreviation A(s±, s 2 , s 3 ) stands for a specific 2-loop tadpole 



N(si,s 2 ,s 3 ) 



2T 2 ((2s 1 23 - 2d) 

P Q [Q 2 ]^[P 2 ] s HiP-Q) 2 ) S3 (2vrT)2 
1 



A( Sl ,S 2 ,S 3 ) = t r ^ r ^„Jr, ^21s 3 = (9„TVsi2 3 -2d N ( S ^ S ^^), (5.5) 



\p 2 + l] s i[q 2 + l] s i[(p - q) 2 ] s 3 



pq 

r(s 13 -i)T(s 23 -l)T^-s 3 )T(s 12 3-d) 

(47r) rf r( Sl )r( S2 )r(d/2)r( Sl2 33 - d) ' 

In Eq. (5.2) we make use of the one-loop subtracted quantities 

P_ 2h 

p) 2 [p 2 y p 2 ' 



(5.6) 



(5.7) 



An(Q) " % [R 2 ) e (R - Q) 2 WF^ 1 WY 1 (5 ' 8) 



as well as the three pieces 



B i AU(P)Q A o I Bii-f'j An ( p )^c 

NM ~%r Q [Q 2 ] N {P-Q) 2, ~ jpjq [Q 2 ] N (P - Q) 2 ' 1 j 

B% M = rf>^l, (5,0) 

where the primed sums denote Yl' n = Xm^cr ^ turns out, however, that BjJ contains an 
additional IR divergence which can be taken into account either by means of IBP reduction 
[15] or by subtraction by hand [12, 16], adding the appropriate zeros (massless tadpoles which 
vanish in dimensional regularization). Performing a transformation to coordinate space in 
d = 3 dimensions and evaluating the remaining sums give the Id integral representations 



B n,m\ €=0 = r(jV) (47r)2jV i^ T ) M I drB^ M (r)An(r), (5.11) 
P>N,M( r ) = E c ^r N ^ {Lbv_ 2+i _ M (e- 2r ) + coth(r) Li N . l+t _ M (e~ 2r )} , (5.12) 



i=0 

N-2N-2+n 



^) = -E E f^^^^'^" 3 ^^-^ 6 "^' (5 - 13) 
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N-2 r 7V— 1 — ' 1 1 

Bn,m(t) = E c ^ r ^ 3 " ~^T^ U N . 1+l .M(e~ 2r ) - - U N . 2+l . M (e- 2r ) , (5.14) 
with A-7r(r) = coth(r) — - — | and where the cjv,j are Fourier coefficients given by 

max(s— 2,1— s) 



— e^/a-Cm) = £ ^ (5.15) 

n=0 

and aAr ;n can be obtained by IBP reduction of the inner sum-integral of Eq. (5.9), see [15]. 

Putting all ingredients together for the special case B^^ (which is needed for ?n|, being 
the coefficient of in Eq. (C.14)), evaluating the finite pieces numerically, 

S&L, = {Li„Ke-0+(coth(r) + i)Lio(e--) + ^Li l( e--) 

x 0.029779678110507967168(1) , (5.16) 



2 (4vr) 4 

T 2 f 00 f 1 1 1 

r 2 

x 0.0020065925001817061293(1) , (5.17) 



2 (4tt) 4 

and using (from IBP, see Eq. (5.4)) 

e ,(3,U.2,0,d) = (rf _ 2) ^-_ 4 | (rf _ 7) 4, 2 . (5.18) 

we obtain as final result for this new master integral (with Z[ = £'(— 1)/£(— 1)) 
T 2 (4ttT 2 )- 3£ 



5 



3,2 



32 (4vr) 4 e 2 



1 + + 7E + 2Z(J e + 70.32026114816592109(1) e 2 + C(e 3 

(5.19) 

For an important cross-check of this result, see App. B. 



6 Applications 

To emphasize the necessity to pursue the matching computations as outlined in this note, let 
us briefly discuss two applications that would become relevant once full results are available. 

The first immediate application involves the Debye screening mass mj| of Sec. 2.1 and 
concerns higher-order perturbative contributions to basic thermodynamic observables, such 
as the pressure of hot QCD. In fact, once the quantity I1e3(0) of Eq. (C.14) has been fully 
determined, the mass term of EQCD (cf. Eq. (2.4)) is available at NNLO, m| ~ g 2 T 2 [l + g 2 + 
g 4 + 0(<? 6 ) ], where g is the dimensionless gauge coupling of full QCD. Now, in the context of 
the effective theory setup for hot QCD, it turns out that the lowest-order EQCD contribution 
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to the full pressure, coming from the quadratic part of £eqcD) enters as ~ Tm| [5], which 
translates to T 4 g 3 [l+g 2 +g 4 + 0(g 6 ) ], such that our 3-loop coefficient contributes to 0(g 7 ) in 
the QCD pressure. According to the systematics of effective theory, due to the fact that there 
are typically large logarithms, a systematic g e evaluation of the pressure (almost completely 
known at present, only missing a well-defined perturbative 4- loop computation [5, 17]) has 
actually been coined physical leading order, since it is the first order where all three physical 
scales (hard/soft/ultra-soft) have contributed. In this respect, the 0(g 7 ) term would simply 
be next-to-leading order, and allow for a first serious investigation of convergence properties. 

Leaving the incomplete 0(g 6 ) (for which there exist numerical estimates, however, from 
comparisons with lattice data, see e.g. [18]) aside for the moment, there are other sources of 
0(<? 7 ) contributions, of course: from the the MQCD pressure plus NLO matching of the 3d 
MQCD gauge coupling g^; from the terms proportional to the quartic coupling Ae in the 
3-loop EQCD pressure; from the 5-loop EQCD pressure (at Ae = 0), which entails one of the 
conceptually simplest (3d, super-renormalizable, massive, vacuum-diagram) computations at 
the 5-loop level, for which techniques are presently developed by several groups; and from the 
leading terms of some higher-order operators in the EQCD Lagrangian, denoted by dots in 
Eq. (2.4), but classified in [19]. All but the last two of these additional g 7 contributions are 
already known. 

A second immediate application involves the 3d EQCD gauge coupling g\ of Sec. 2.2 
and concerns precision-tests of the dimensional reduction setup, such as for the spatial string 
tension a s , which parameterizes the large-area behavior of rectangular spatial Wilson loops. 
As has been demonstrated in Ref. [7], it can be systematically determined, as a function of 
the temperature T, in the dimensionally reduced effective theory setup, and then compared 
to non-perturbative 4d lattice measurements. It turned out that the NLO result for gj| as 
obtained in [7] represents a considerable improvement over a 1-loop comparison - giving a 
sizable correction factor as well as a first estimate of (renormalization) scale dependence - 
while leaving room for NNLO effects, for which our 3-loop result for il^ 3 (0) of Eq. (C.15) is 
the last missing building block. 

7 Conclusions 

We have successfully reduced the NNLO contributions to the matching parameters m| and 
g\ to a sum of scalar sum-integrals. These matching parameters play an important role 
in higher-order evaluations of basic thermodynamic observables and in precision-tests of the 
dimensional reduction setup respectively, and hence are needed with high accuracy. Our result 
passes the non-trivial checks of transversality as well as gauge-parameter independence. 

In a next step, a number of master integrals have to be evaluated. Although we managed 
to map all of them to the relatively simple class of basketball-type ones, the somewhat large 
number of masters which we need demand a semi-automated evaluation strategy, which still 
has to be developed. As a first and encouraging step towards this goal, we have demonstrated 
a systematic method to evaluate a certain class of such basketball-type sum-integrals. 
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Once full results for the matching coefficients discussed here become available, there are 
immediate applications to quantities of phenomenological interest, such as the pressure of hot 
QCD, or the spatial string tension, as discussed in Sec. 6 above. However, these concrete 
applications will have to await progress in the art of sum-integration for now. 



Acknowledgments 

This work was supported by the Deutsche Forschungsgemeinschaft (DFG) under contract 
no. SCHR 993/2, by the BMBF under project no. 06BI9002, and by the Heisenberg Pro- 
gramme of the DFG. contract no. SCHR 993/1. 



A One- and two-loop vacuum sum-integrals 

The one-loop bosonic tadpole is known analytically and reads 

2^ 3 / 2 T 4 / /x 2 \ e r(m-| + e) 



((2m - n - 3 + 2e) , (A.l) 



r =4 p °" _ 

m ~ Tp (P 2 ) m (2^T) 2m "« \ttT 2 J r(m) 
whereas the fermionic tadpole can be related to the corresponding bosonic one via 



f pn 

jri — /<)2m- n— 3+2e i \ jn l \ n\ 

As mentioned above, via integration-by-parts relations all two-loop integrals are expressible 
in terms of products of two one-loop tadpoles which means they are also available analytically 
up to arbitrary order in e. 



B Check of new sum-integrals 



We can cross-check our new result given in Sec. 5 using IBP reduction of the V-type topology 
which gives 



V = I 



0,0,0 



4-E?3,2 -I — -r £?2 ,0 



'1,1,1,1,1,0; 0,0,0" 3(d _ 3) 2 

where V stands for the spectacles-type diagram given in [13]: 

1 



2(d-4) 



(B.l) 



V = 



PQR F*Q*(P-Qy&(P-R)* 
T 2 (47rT 2 e 7E )- 3e 



4(4vr) 4 e 2 



\ + AK 1 e + AK 2 e 2 + 0{e^) 



(B.2) 



with AK\ = | + 47e + 2Z{, while AK 2 is known only numerically. Writing the coefficients of 
our basketball-results, given in Eq. (5.19) above as well as Eq. (26) of [15], as 



2\-3e 



B. 



3,2 



T 2 (4^T 2 ) 
32(4vr) 4 e 2 



3320 



+ 632ie + &322e 2 + C(e 3 ) 



(B.3) 
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B 



T 2 {4ttT 2 ) 



l + b 2 ie + b 22 e 2 + 0(e 3 ) 



2,0 8(4vr) 4 e 2 

to match the leading term of V it follows that the linear relations 



^320 



1 



^321 



621 + 4 



^322 



b 22 + 46; 



21 



(B.4) 



(B.5) 



have to be satisfied. Our results presented above do indeed confirm these relations, which we 
take as a nice check of our generic parameterizations. Eq. (B.5) provides a welcome check of 
our numerical constants. 



C Expansion coefficients up to three loops 

For convenience, we here repeat the one- and two-loop coefficients that were already computed 
in [7], adding the second derivatives that are needed for Eq. (2.8). The one-loop coefficients 
up to second derivative read 



n T i(o) = 0, 

n E i(o) = (d - 1) [c A (d - 1)/? - 2n41 



L 2 i 



n^(o) 



C A 



41 1 
10 ~ 10 



d + 2£- V 



15 



3 - 



n£i(o) = % 



0"i 7 1 < 2 



The two- loop coefficients up to first derivative are given by (see also [7]) 

n T2 (o) = 0, 

n E2 (o) = (d- i)(d -3)((i + o \mii -(d- i)c A i?l c A / 2 °+ 



n' T2 (o) 



+ 2 N { C F 
(d-3)(d-4) 



fo 



2 r 1 



(d - 7)(d - 5)(d - 2)d 



(-14-42d + 8cf )C A / 2 °/ 2 ° 



4C F + (1 - 6d + d 2 )C A JVfiJ ^2 



y-6d 2 + |d- 



■14 + 41d - 12d 2 + d 3 )C F 



(C.l) 
(C.2) 

(C.3) 
(C.4) 
(C.5) 
(C.6) 

(C.7) 
(C.8) 
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(d-1) 



(144 - 31d + d 2 ) (1 - d)C A Ji + 2iV f l 



3d{d - 7) 
4(d - 6)(d - l)C F iV f [/? - /? 
(d-3) 



4%, 



(C.9) 



n - W = 2(d-7)(d-5)(d-2)d j ^ + 31M " ^ + * 7 * ~ ^ ° l /2 ° /2 ° + 
+ 2(d - 4)(d - 1) [ (2 - 5d + d 2 ) C A + 8C F j N { 1$ I$+ 

+ (d - 1) (24 - 7d 2 + d 3 ) C A - 2 (28 + 2d - 7d 2 + d 3 ) C F N { |g 1$ \+ 
(d-3)£ 



+ 



24(d - 2) 
(d-3)(d-l) 



3 (16 - 13d + 3d 2 ) £ - 4 (44 - 29d + 7d 2 - d 3 ) 
(d-1) 



r*2 to to 

°A J 2 1 2~ 



ec A iVf/ 2 i 2 + 



4(6 + 15d - 10d 2 + d 3 ) 



x C F N f 



7O 



6(d-7)d 

2 (-72 + 42d - 13d 2 + d 3 ) + 2(d - 7)d 2 £+ 



+ (d - 7)(d - 6)d£ 2 l [(d - l)C A /i° - 2^/° 



(CIO) 



For presenting the outcome of the reduction procedure for the three- loop contributions, 
which constitutes the main result of this paper, we denote the master integrals as in Fig. (2), 
i.e. /,/ for the 1-loop tadpoles of Eqs. (A.l) and (A. 2), and 



a,b,c,d,e,f a,b,c,d,e,f;0,0,0 



K 



I 



a,b,c,d,e,f a,6,c,e£,e,/; 0,0,1 



I 



a, ^,7 



- a,p,7 

" ' a,b,c,d,e,f J a,6,c,<i,e,/; 1,1,0 



(C.ll) 
(C.12) 
(C.13) 



are 3- loop basketball- type integrals in a slightly more compact notation than Eq. (3.1). The 
results needed for Eqs. (2.8) and (2.12) then read 



rr m\ /~<3 \ 7OAO 1 t0,0,2 , T 0,2,0 . T 2,0,0 . T l,3,0 

llE3(UJ - C A [«1 J 2 ,1,0,0,1,1 + "2^2,2,0,0,1,1 + a 3J 3 ,l,0,0,l,l + a 4«V 3i i i0 ,0,l,l + a 5 -4,1,0,0,1,1 + 

76,0,0 . 76,4,0 , 77,3,0 . r 70 70 1 70 70 70^ 

a 6 J 5,l,0,0,l,l + "7^5,3,0,0,1,1 + a 8^6,2,0,0,l,l + hh + a ^hhh 

1 ri2 at I 7^0,0,0 . 7^2,0,0 . 7v-l,l,0 . 7^2,0,0 

+G A 7V f [aiiA ljlA021 + ai2ii 110A22 + ai 3 A 1)1)0i0)3jl + auK 1X0fi3>1 + 
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7^7.0,3 . 7^9,0,1 . 7^10,0,0 . 7/7,3,0 

"39-^6,1,0,0,2,1 + "40-^6,1,0,0,2,1 + "41 -^6,1,0,0,2,1 + "42-^6,2,0,0,1,1 + 

7^8,2,0 . 7^8,0,2 i 7^8,2,0 i J^ 9 ,0,1 

"43-^6,2,0,0,1,1 + a 44-K 7) i i0 ,0,l,l + "45-^7,1,0,0,1,1 + "46-^7,1,0,0,1,1 + 

7^9,1,0 . 7^10,0,0 . T 0,0,0 . r 2,0,0 

"47-^7,1,0,0,1,1 + "48-^7,1,0,0,1,1 + "49-^2,1,0,0,1,1 + "50-^3,1,0,0,1,1 + 

asJinil + a^llll + a^Ipl + a^PJPl + a 55 I°P 1 li+ 
«56/i /i°/ 3 + «57/i°/ 2 °/ 2 + a 5 8/i /i /3° + «59/i / 2 °/ 2 + + 

\m at%\ r 0,0,0 . r 0,0,2 1 r 0,2,0 1 r 2 A° 

+C A iV f ["61-^2,1,0,0,1,1 + "62-^2,2,0,0,1,1 + "63^3,1,0,0,1,1 + "64-^ 3 ,i, ,0,l,l + 

T 1,3,0 . r 6,0,0 . r 6,4,0 . r 7,3,0 

"65-^4,1,0,0,1,1 + "66-^5,1,0,0,1,1 + "67-^5,3,0,0,1,1 + "68-^6,2,0,0,1,1 + 



0,0,0 . T 0,0,2 T 0,2,0 . T 2,0,0 

'2,1,0,0,1,1 + "74-^2,2,0,0,1,1 + "75-^3,1,0,0,1,1 + "76-^3,1,0,0,1, 
6.0,0 , fO fO fO , ^„ fO tO tO , fO 7O fol , 



+ 



"77i5;?A0,l,l + "78A°A°4° + "79A°4°4 + "8oA°^ j 

7^0,0,0 . 7^2,0,0 . t-1,1,0 . 7^0,0.0 

["81^1,1,0,0,2,1 + "82-^1,1,0,0,2,2 + a S3 K l,l,0fi,3,l + "84-^2,1,0,0,1,1 + 

7^0,2,0 . 7^2,0,0 . 7,-1,1,0 . 7^2,0,0 

"85-^2,1,0,0,2,1 + "86-^2,1,0,0,2,1 + "87-^2,2,0,0,1,1 + "88-^2,2,0,0,1,1 

7^0,2,0 . 7^1,1,0 , 7^2,0,0 , 7^1,3,0 

"89-^3,1,0,0,1,1 + "90-^3,1,0,0,1,1 + "91^3,1,0,0,1,1 + "92-^3,1,0,0,2,1 

7^0,4,0 . 7^1,3,0 . 7^2,2,0 , 7^4,0,0 

"93-^-3,2.0,0,1,1 + "94-^-/1 1 n n 1 1 + "95-^-/1 1 n n 1 1 + "9fi-"^ 1 n n 1 1 



1,3,0 , 7>-2,2,0 . 7^4,0,0 

4,1,0,0,1,1 + "95- n -4,l,0,0,l,l + "96-^-4,1,0,0,1,1" 

3,3,0 . 7^5,1,0 . r 0,0,0 

5,1,0,0,1,1 + "99 A 5,1,0,0,1,1 + "100-^2,1,0,0,1,1 

rO fO fO , _ _ tO tO T , _ „ . fO fO tO, 



+ 
+ 
+ 



+C A N { C F 



7^6,0,0 . 7^3 

"97^4,2,0,0,1,1 + "98^5,,, 
"101-^3+0,0,1,1 

+ ai 02 /i U /i U / 3 + "103-^1/3 + "104/1 / 2 U / 2 U + 

ai05/i°/2°/2 + "l06/l 2 /2 / 3 ° + "l07/l°/l°/3 + "lOs/?^ + a m I^^+ 

ano/i /2 / 2 + "iii/i 2 / 2 °/3] + 

7^0,0,0 . 7^2,0.0 . 7^0.0,0 . 7^0,2,0 

"112-^1,1,0,0,2,1 + "113-^1,1,0,0,2,2 + "114-^2,1,0,0,1,1 + "115^2,1,0,0,2,1 + 

7^2,0,0 . 7^1,1,0 . 7v- 2 ,0,0 . 7^0,2,0 

"116-^2,1,0,0,2,1 + "117-^2,2,0,0,1,1 + "118^2,2,0,0,1,1 + "H9 K 3,1,0,0,1,1 + 

7^1,1,0 . 7^2,0,0 . 7^1,3,0 . 7^0,4,0 

"120-^3,1,0,0,1,1 + "121-^3,1,0,0,1,1 + "122-^3,1,0,0,2,1 + "123-^3,2,0,0,1,1 + 

-1,3,0 . 7^2,2.0 . 7^4.0,0 . 7^6,0,0 

a 1 n n 1 1 + a-\9F,J\ A 1 n n 1 1 + ttl^fiii-zi 1 n n 1 1 + Ctl97 1\ a o n n i i + 



,U,U,1,1 0,X,U,U,X,1 0,-L,U,U,Z,l 1 

0124^4,1,0,0,1,1 + "125-^4,1,0,0,1,1 + "126-^4,1,0,0,1,1 + "127-^4,2,0,0,1,] 

7^3,3,0 . r 0,0,0 . 7 2,0,0 . fOfOrO, 

"128-^5,1,0,0,1,1 + "129-^2,1,0,0,1,1 + "130-^3,1,0,0,1,1 + a m I l 1 lh + 

"l32/l /2 / 2 ° + "l33/l / 2 °/2 + "l34/l 2 / 2 °/3 + "iss/?^^ + "136^?^ 
72 fO 70" 



n^ 3 (o) 



"l37/l 2 / 2 °/3 

' o t-0,0,0 . q t0,0,0 , a 70,0,2 , R t0,2,0 , 

Pi J 2, 2,0,0, 1,1 + P2 J 3,l,o,0,l,l + ^3^3,2,0,0,1,1 + "4^4,1,0,0,1,1 + 

o t2,2,0 , a t-4,0,0 R 7 8,0,0 , R 7 7,3,0 

^5^5,1,0,0,1,1 + P6^5, 1,0,0, 1,1 ' P7*^7,l,0,0,l,l + P8<^7,3,0, -1,1,1 + 



(C.14) 



A/; 



8,2,0 
8,2,0,- 



ua + /3io/i /i / 4 + Ai/i U / 2 U / 3 U + /5i2/ 2 u / 2 U /2 U 



7O 7O 



70 70 70 



+ 



^i, a 7^0,0,0 . a r^0,0,0 . o 7^2,0,0 . a T A,l,0 

+0 A iV f Pl3-Ki,i,o,0,2,2 + P14-K 1,1,0,0,3,1 + ^15-^1,1,0,0,3,2 + Pl6- K l,l,0,0,4,l + 
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o 7,-2,0,0 , a 7^0,0,0 , a 7^0,2,0 , n 7^2,0,0 

Pl7 K l, 1,0,0,4,1 + ^18-^2,1,0,0,2,1 + ^19-^2,1,0,0,3,1 + P20-K 2 ,l,0,0,3,l 



a ^4,0,0 o 7^0,0,0 o ^0,0,0 fl ^0,2,0 

^21-^2,1,0,0,4,1 + P22-K 2 ,2,0,0,l,l + P23-K 3 ,l,0,0,l,l + P24-K 3 ,i, ,0,2,l 

ion Ann r\ n o 



+ 
+ 

^25-^3,1,0,0,2,1 + P26-K 3 ,l,0,0,3,l + P27-K 3 ,l,o,0,3,l + P28-K 3 ,2,0,0,l,l + 

1 + 
i + 



^37^2,0,0, 



* ? } ! ) 15 7 51 5 5 5 5 5 

o t>-0,2,0 . o 7^1,1,0 , o 7^0,0,2 . o r^0,2,0 

P29 A 3,2,0,0,1,1 + ^30^3,2,0,0,1,1 + P31- K 4,l,0,0,l,l + ^32^4,1,0,0,1,1" 

'° _l_ R r^2,0,0 . n 7^1,3,0 . a 7^0,4,0 

,0,0,1,1 + P34-"-4 i i j o,0,l,l P35-"- 4,1,0,0,2,1 ~+~ P36 A 4 I 2,0,0,1,1 
',0 . a 7^1,1,2 . o tv-1,3,0 . a v-2,0,2 

,0,0,1,1 + P38-tt 5 ,l,0,0,l,l + P39-K5,l,0,0,l,l + P40-K5,l,0,0,l,l - ' 

>° _L R IT 3 ' 1 . i R 1^4,0,0 i R 1^3,3,0 , 

'"^5,1,0,0,1,1 + P43-K 5 ,l,0,0,l,l + P44-K 6! i )0 ,0,l,l+ 

^ 6 '°.° , . , a ^8,0,0 _ _ _ , ~5,3,0 _ _ , 



+ 



1,1 + 
+ 



-C A Nf 



+NfC ¥ 



+N f C% 



P41-^5,'l,0,0,l,l + P42— 

P45- K 6,l,0,0,l,l + P46- K 6,l,0,0,l,l + P47-K 6 ,2,0,0,l,l + P48- K 7,l,0,0,l,l + 
a z>-6,0,2 . a 7^7,1,0 . o 7^8,0,0 . o 7^8,2,0 

P49- K 7,l,0,0,l,l + P50-«-7,l,0,0,l,l + P51 A 7,1,0,0,1,1 + P52- K 7,2,0,0,l,l + 

^53^7^2,0,0,1,1 + P54-^8,'l,'-l,0,2,l + ^55-^8,1,-1,0,2,1 + /^56-^8,'l,0,0,l,l 
a 7^9,0,1 . a 7^9,1,0 . n 70-10,0,0 . n jo-8,1,1 

P57-K 8 ,l,0,0,l,l + P58-K 8 ,l,0,0,l,l + P59-K 8 ,l,0,0,l,l + £60^8,2,-1,0,1, 

_L R r 0.0.0 jo r0,0,0 , o r 4,0,0 

,1 + P62 ^2,2,0,0,1,1 + P63i'3 ) i ) o l O,l,l + P64%,1,0,0,1,1 

p^nnii + pw?/^ + + a»mj§ + 

+ /Jn^/o/o + p 72 p 2 m + ^i?/?^ + /3 74 / 1 o /T/4 o + 

/W?/^ + + /?77/ 1 °/l /4° + + 

~ a r°> ' i a fOflfi . n r 0,0,2 n r 0,2,0 

P84-^2,2,0,0,l,l + P85-^3,l,0,0,l,l + P86-^3,2,0,0,l,l + P87-^4,l,0,0,l,l + 

o r 2,2,0 n r 4,0,0 , o r 8,0,0 n r7,3,0 

^88-^5,1,0,0,1,1 + ^89^5,1,0,0,1,1 + ^90-^7,1,0,0,1,1 + P91 ^7,3,0,-1,1,1+ 

^92-^8^0,-1,1,1 + ^93-^1-^1 ^4 + PV^l I1I4 + /?95-^l-^2-^3 + P96-^l ^2 ^3 + 
^97-^1-^2-^3 + PmHH^ + /?99-^2-^2-^2 + 

/3l00-C'2i2A0,l,l + ^lOl^^O.O.l,! + / 3 102-^3!2A0,1,1 + / 5 103^1AO,l,l + 

fl r 2,2,0 , o r 4,0,0 , o r 8,0,0 , fl fO fO fO , 

Pl04-^5,1,0,0,1,1 + P105-^5,1,0,0,1,1 + Pl06-^ 7 ,1,0,0,1,1 + + 

^108-^1-^2-^3 + /3l09-^l-^2-^3 + AlO-^l-^2-^3 + ^111-^2-^2^2 + 

"o tv-0,0,0 . o 7^0,0,0 . o 7^1,1,0 . o k-2,0,0 

^112-^1,1,0,0,2,2 + PII3A1, 1,0,0,3,1 + PH4-Kl, 1,0,0,4,1 + ^115^1,1,0,0,4,1 + 

o 7>-0,0,0 . o 7^0,2,0 . a 7^2,0,0 . „ r -0,0,0 

Pll6-K 2 ,i, ,0,2,l + Pll7-K 2 ,i, ,0,3,l + PH8- K 2,1,0,0,3,1 + Pll9 A 2,2,0,0,l,l + 

o 7>-0,0,0 , o ^0,2,0 , n 7^1,1,0 , n 7^1,3,0 , 

Pl20-ft- 3 ,l,o,0,l,l + ^121-^3,1,0,0,2,1 + ^122^3,1,0,0,2,1 + Pl23-K 3 ,l,0,0,3,l + 



/3l24<.„ 

o t-1,1,0 , a t-2,0,0 , o 7>-l,3,0 , a tv-0,4,0 

Pl28-K 4;1)0 ,o,l,l + Pl29-K 4 ,l,0,0,l,l + Pl 3 A 4,l,0,0,2,l + Pl31- K 4,2,0,0,l,l 

o z^-4,0,0 . o 7^1,3,0 . n 7^2,2,0 . „ tv-3,1,0 

Pl32-«- 4 , 2 ,o,0,l,l + Pl33-K 5 ,l,0,0,l,l + P134-K 5 ,1,0,0,1,1 + Pl35-K 5 ,l,0,0,l,l 

o t>-4,0,0 . o 7^3,3,0 . n k-4,2,0 . „ 7^6,0,0 

Pl36-K 5 ,l,0,0,l,l + Pl37-K 6 ,l,0,0,l,l + Pl38-K 6 ,l,0,0,l,l + P139 A 6,1,0,0,1,1 



+ 
+ 
+ 
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/Wl^ + /Wl"^ + /3l50/ll2 /3 + PisJIlP! + /3l5 2 /l 2 /3 4° + 
/3l53-f 2 -f 2 -^ 2 + ^154-^1^1-^4 + ^155-^1 -^2 ^3 + ^156-^1 ^2 ^3 + ^157-^1^1-^4 + 

/3l58-^l-^ 2 -^3 + A 59 ^1-^4 + P160I1I3I3 + ^161-^2^2-^2 + ^162^2 ^2 ^2 + 

+C A A r fC , F /^leS^X'o.O^^ + / 3 164-^°f j o,0,3,l + ^leS-^f^o.CS^ + ^166^1,1,0,0,4,1 + 

o tv-0,0,0 , o 7^0,2,0 , q 7^2,0,0 , o t^O.0,0 , 

^167-^2,1,0,0,2,1 + Pl68-K 2 ,l,0,0,3,l + ^169^2,1,0,0,3,1 + W0-K 2 ,2,0,0,l,l + 

n n n non iin 1 o n 



-<3 -r Pl98- l l-'2- t 3 P199 J 1 J 2 J 3 



/3l95-^5j'o,0,l,l + ^196-^1^1/4 + #197-7? ^2 J _ 

#200 -^1-^4 + P201I1I3I3 + #202 -^2 -^2 -^2 + #203-^1 -^1^4 + #204-^1 ^2 ^3 + 
#205-^1-^2^3 + ^206-^1-^2-^3 + #207-^1 ^2 ^4 + #208^1 ^3-^3 + #209-f 2 ^2 ^2 + 

/3 2 io/ 2 / 2 4 ] • (C15) 

Looking at the master integrals that are needed for the above two lengthy expressions let 
us note that, while most of them have factors of -Po etc- i n the numerator, only eight of them 
(those multiplying #{8,9,54,55,60,61,91,92}) contain irreducible scalar products in the numerator 
and hence need methods for their evaluation that go beyond those presented in Appendix 5 
(see, however, Ref. [11, 20], where examples of such sum-integrals were treated). Also, some 
of the masters (such as e.g. those multiplying a{7,s}i #{8,9}) involve somewhat large powers 
of propagators, which is a consequence of our ordering prescription. However, as was shown 
in App. 5 in terms of the generic power N, this does not seem to be a particularly difficult 
obstacle. 

We refrain from listing the coefficients ai...i37 an d #i...2io here. They have the general 
form X^ n £™ Pn(d)/q n (d), where £ is the gauge parameter (see Eq. (2.14)) and p, q are polyno- 
mials in d. The full expressions for Eqs. (C.14) and (C.15) are provided in computer-readable 
form on [21]. 
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